Chapter 4 

Fluid Description of Plasma 



The single particle approach gets to be horribly complicated, as we have seen. 

Basically we need a more statistical approach because we can't follow each particle separately. 
If the details of the distribution function in velocity space are important we have to stay 
with the Boltzmann equation. It is a kind of particle conservation equation. 



4.1 Particle Conservation (In 3-d Space) 



Av 




Ax 

Figure 4.1: Elementary volume for particle conservation 

Number of particles in box AxAyAz is the volume, AV = AxAyAz, times the density n. 
Rate of change of number is is equal to the number flowing across the boundary per unit 
time, the flux. (In absence of sources.) 

d 

— —-[AxAyAz n] = Flow Out across boundary. (4.1) 

Take particle velocity to be v(r) [no random velocity, only flow] and origin at the center of 
the box refer to flux density as n\ = J. 

Flow Out = [J,{0,0,Az/2) - J,{0,0,-Az/2)] Ax Ay + x + y . (4.2) 
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Expand as Taylor series 

U0,0,V) = UO) + ^Jz.V (4.3) 

oz 



So, 



d 

flow out ^ —(nVz)AzAxAy + x + y (4.4) 
oz 



= AW . (nv). 

Hence Particle Conservation 



d 

—n = V.(nv) (4.5) 



Notice we have essential proved an elementary form of Gauss's theorem 



V.Ad'r = / A.dS. (4.6) 

Jd-y 

The expression: 'Fluid Description' refers to any simplifled plasma treatment which does 
not keep track of v-dependence of / detail. 

1. Fluid Descriptions are essentially 3-d (r). 

2. Deal with quantities averaged over velocity space (e.g. density, mean velocity, ...). 

3. Omit some important physical processes (but describe others). 

4. Provide tractable approaches to many problems. 

5. Will occupy most of the rest of my lectures. 

Fluid Equations can be derived mathematically by taking moment^ of the Boltzmann Equa- 
tion. 

0*'' moment / d^v (4.7) 
1st moment / vrf^f (4.8) 
2nd moment / vvd^f (4.9) 

These lead, respectively, to (0) Particle (1) Momentum (2) Energy conservation equations. 
We shall adopt a more direct 'physical' approach. 
^They are therefore sometimes called 'Moment Equations.' 
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4.2 Fluid Motion 



The motion of a fluid is described by a vector velocity field v(r) (which is the mean velocity 
of all the individual particles which make up the fluid at r). Also the particle density 
n{r) is required. We are here discussing the motion of fluid of a single type of particle of 
mass/charge, m/q so the charge and mass density are qn and mn respectively. 

The particle conservation equation we already know. It is also sometimes called the 'Conti- 
nuity Equation' 

|-n + V.(nv) = (4.10) 

at 

It is also possible to expand the V. to get: 

d 

— n + (v.V)n + nV.v = (4.11) 

ot 

The signiflcance, here, is that the first two terms are the 'convective derivative" of n 

D d d , , 



so the continuity equation can be written 

D 



n = -nV.v (4.13) 
Dt ^ ' 



4.2.1 Lagrangian & Eulerian Viewpoints 

There are essentially 2 views. 

1. Lagrangian . Sit on a fluid element and move with it as fluid moves. 




Figure 4.2: Lagrangean Viewpoint 



2. Eulerian . Sit at a flxed point in space and watch fluid move through your volume 
element: "identity" of fluid in volume continually changing 

^ means rate of change at fixed point (Euler). 

^ = J: = ^ + v.V means rate of change at moving point (Lagrange). 
v.V = |!|- + ^|- + l!|-: change due to motion. 
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Figure 4.3: Eulerian Viewpoint 



Our derivation of continuity was Eulerian. From the Lagrangian view 



D _ d AN __AN d _ 1 dAV 

Di^~ dt AV ~ ~AV^ Jt ~~^AV~dr 



(4.14) 



since total number of particles in volume element (AA^) is constant (we are moving with 
them). (AV ^ AxAyAz.) 



Now —AV 
dt 



dAx , , dAy , , dAz , , 
AyAz + -^AzAx + -^AyAx 



dt 
AV 



dt 



dt 



1 dAx 1 dAy 1 dAz 

+ H 

Ax dt Ay dt Ax dt 



(4.15) 
(4.16) 



Hence 



and so 



But 



d{Ax) 
dt 



{Ax/2) - {-Ax 12) 

...y ...z 



Ax -p— etc. 
ox 



dt 



AV = AV 



dvx ^ dvy ^ dvz 
dx dy dz 

D 



AVV.Y 



Dt 



n — — nV.v 



(4.17) 
(4.18) 

(4.19) 
(4.20) 



Lagrangian Continuity. Naturally, this is the same equation as Eulerian when one puts 



+ v.V. 



D_ — d_ 
Dt dt 

The quantity — V.v is the rate of (Volume) compression of element. 



4.2.2 Momentum (Conservation) Equation 

Each of the particles is acted on by the Lorentz force g[E + Ui AB] (u^ is individual particle's 
velocity). 

Hence total force on the fluid element due to E-M fields is 



(g [E + Ui A B]) = AA^ 5 (E + V A B) 



(4.21) 
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(Using mean: v — u/AA?".) 

E-M Force density (per unit volume) is: 

Fem =ng(E + v AB). (4.22) 

The total momentum of the element is 

J2 = m AN V = AV mnv (4.23) 

i 

SO Momentum Density is mnv. 

If no other forces are acting then clearly the equation of motion requires us to set the time 
derivative of mn\ equal to Fem- Because we want to retain the identity of the particles 
under consideration we want D/Dt i.e. the convective derivative (Lagrangian picture). 

In general there are additional forces acting. 

(1) Pressure (2) CoUisional Friction. 



4.2.3 Pressure Force 

In a gas p{— nT) is the force per unit area arising from thermal motions. The surrounding 
fluid exerts this force on the element: 



Figure 4.4: Pressure forces on opposite faces of element. 



Net force in x direction is 

- p{Ax/2)AyAz+p{-Ax/2)AyAz (4.24) 

~ -AxAyAz ^ = -AV ^ = -AV (Vp), (4.25) 
ox ox 

So (isotropic) pressure /orce density (/unit vol) 

Fp = -Vp (4.26) 

How does this arise in our picture above? 

Answer : Exchange of momentum by particle thermal motion across the element boundary. 

Although in Lagrangian picture we move with the element (as defined by mean velocity v) 
individual particles also have thermal velocity so that the additional velocity they have is 

Wj = Uj — V 'peculiar' velocity (4.27) 
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Because of this, some cross the element boundary and exchange momentum with outside. 

(Even though there is no net change of number of particles in element.) Rate of exchange 

of momentum due to particles with peculiar velocity w, (Pw across a surface element ds is 

/(w)mw d^w X w^_^ds (4.28) 
mom"* density at w rate across ds 
Integrate over distrib function to obtain the total momentum exchange rate: 

ds. j mww/(w)(i^w (4.29) 
The thing in the integral is a tensor. Write 

p = J mwwf{w)d^w (Pressure Tensor) (4.30) 
Then momentum exchange rate is 

p . ds (4.31) 
Actually, if /(w) is isotropic (e.g. Maxwellian) then 

Pxy — J m Wx Wy f{w)d^w — etc. (4.32) 

and = j mwl f{w)d^w = nT{= pyy = p^^ = 'p') (4.33) 
So then the exchange rate is pds. (Scalar Pressure). 

Integrate ds over the whole AV^ then x component of mom"* exchange rate is 

'Ax\ . . /—Ax^ 

and so 

Total momentum loss rate due to exchange across the boundary per unit volume is 

Vp (= -F,) (4.35) 

In terms of the momentum equation, either we put Vp on the momentum derivative side or 
Fp on force side. The result is the same. 

Ignoring Collisions, Momentum Equation is 

(mnAVv) = [Fem + Fp] AV (4.36) 
Recall that nAy = AA^ ; g(AA^) = 0;so 

Dv 

L.H.S. = mnAV — . (4.37) 
dt 



p 



^) AyAz - -p AyAz = AV (Vp), (4.34) 



Thus, substituting for F's: 

Momentum Equation. 



Dv dw \ 

mn— = mn — + v.Vv = on (E + v A B) - V» (4.38) 
Dt \ at 
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4.2.4 Momentum Equation: Eulerian Viewpoint 

Fixed element in space. Plasma flows through it. 

1. E.M. force on element (per unit vol.) 

Fem = nq{E + w AB) as before. (4.39) 

2. Momentum flux across boundary (per unit vol) 

= V. J m{v + w){v + w) f{w) d^w (4.40) 

= V.{jm{vv+ y w + wv +ww) /(w) li^w} (4-41) 
integrates to 

= V.{mnvv + p} (4.42) 

= mn(v.V)w + mv [V.(nv)] + Vp (4.43) 

(Take isotropic p.) 

3. Rate of change of momentum within element (per unit vol) 

d 

= —{mnv) (4.44) 

Hence, total momentum balance: 

d 

— {mnv) + mn(v.V)v + mv [V.(nv)] + Vp = F^m (4.45) 
Use the continuity equation: 

^ + V.(nv) = , (4.46) 

to cancel the third term and part of the 1st: 

d dn dv dv 

— {mnv) + mv (V. {nv)) — mv{ — + V. {nv)} + mn— = mn— (4.47) 

Then take Vp to RHS to get final form: 
Momentum Equation: 



mn 



dv 

g^ + (v.V)v 



ng (E + V A B) - Vp . (4.48) 
As before, via Lagrangian formulation. (Collisions have been ignored.) 
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4.2.5 Effect of Collisions 



First notice that like particle collisions do not change the total momentum (which is averaged 
over all particles of that species). 

Collisions between unlike particles do exchange momentum between the species. Therefore 
once we realize that any quasi-neutral plasma consists of at least two different species (elec- 
trons and ions) and hence two different interpenetrating fluids we may need to account for 
another momentum loss (gain) term. 

The rate of momentum density loss by species 1 colliding with species 2 is: 

i/i2nimi(vi - V2) (4.49) 
Hence we can immediately generalize the momentum equation to 
(9vi 



+ (vi.V) vi 



dt 

With similar equation for species 2. 



riiqi (E + vi A B) - Vpi - Ui2nimi (vi - V2) (4.50) 



4.3 The Key Question for Momentum Equation: 

What do wc take for pi 

Basically p = nT is determined by energy balance, which will tell how T varies. We could 
write an energy equation in the same way as momentum. However, this would then contain 
a term for heat flux, which would be unknown. In general, the A;*'* moment equation contains 
a term which is a (A; -|- 1)*'* moment. 

Continuity, 0*'* equation contains v determined by 
Momentum, 1** equation contains p determined by 
Energy, 2^'^ equation contains Q determined by . . . 

In order to get a sensible result we have to truncate this hierarchy. Do this by some sort of 
assumption about the heat flux. This will lead to an 
Equation of State: 

pn~'^ — const. (4-51) 

The value of 7 to be taken depends on the heat flux assumption and on the isotropy (or 
otherwise) of the energy distribution. 

Examples 

1. Isothermal: T = const.: 7 = 1. 

2. Adiabatic/Isotropic: 3 degrees of freedom 7 = §■ 

3. Adiabatic/1 degree of freedom 7 = 3. 
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4. Adiabatic/2 degrees of freedom 7 = 2. 



In general, n{i/2)5T = —p{5V/V) (Adiabatic £ degrees) 



I 5T -6V 6n 

So 

p n T \ tj n 

i.e. 

const. (4.54) 

In a normal gas, which 'holds together' by collisions, energy is rapidly shared between 3 
space-degrees of freedom. Plasmas are often rather coUisionless so compression in 1 dimension 
often stays confined to 1-degree of freedom. Sometimes heat transport is so rapid that the 
isothermal approach is valid. It depends on the exact situation; so let's leave 7 undefined 
for now. 



4.4 Summary of Two-Fluid Equations 

Species j 
Plasma Response 



1. Continuity: 



dt 

2. Momentum: 



^''^■+V.Kv,) = (4.55) 



a/ + (v,.v)v, 



= njQj (E + Vj A B) — Vpj — Vj^rij-rrij {vj — v^) (4.56) 



3. Energy /Equation of State: 

PjuJ'^ — const.. (4-57) 

(j = electrons, ions). 



Maxwell's Equations 



V.B = V.E = p/eo (4.58) 
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With 



j = ge^eVe + qiUiWi = 6 {^71^^^ + ZUiVi 

— —eue (vg — Vjr) (Quasineutral) 



(4.60) 
(4.61) 
(4.62) 



Accounting 



Unknowns 

ne,ni 2 



Pe,Pi 

E,B 



6 
2 
6 
16 



Equations 
Continuity e, i 
Momentum e, i 
State e, i 
Maxwell 



2 
6 
2 
8 
18 



&ut 2 of Maxwell (V. equs) are redundant because can be deduced from others: e.g. 

d 



V.(VAE) = 
and V. (V A B) = 



dt 



(V.B) 



1 d 

/^oV.j + ^ (V.E) 



1 ^ 
^ di 



Co 



+ V.E 



(4.63) 
(4.64) 



dt ' ' c 
So 16 equs for 16 unknowns. 

Equations still very difficult and complicated mostly because it is Nonlinear 

In some cases can get a tractable problem by 'linearizing'. That means, take some known 
equilibrium solution and suppose the deviation (perturbation) from it is small so we can 
retain only the 1st linear terms and not the others. 



4.5 Two-Fluid Equilibrium: Diamagnetic Current 

Slab: i ^ i = 0. 

Straight B-field: B = Bz. 
Equilibrium: | = {E ^ -V0) 

CoUisionless: — > 0. 
Momentum Equation(s): 

^i^i(vi-V)vj = njgj(E + A B) - Vpj (4.65) 

Drop j suffix for now. Then take x, y components: 

d / \ dp 

mnv^—v^ = nq{Ex + VyB)-— (4.66) 

mn Vx^Vy — nq{0 — v^B) (4-67) 

LLJb 
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Eg |4.67| is satisfied by taking = 0. Then [466 



I.e. 



or, in vector form: 



nq{E, + VyB) - ^ = 0. (4.68) 

ax 



EAB Vp B 

nq B^ 

EAB drift Diamagnetic Drift 



Notice: 



• In magnetic field (_L) fluid velocity is determined by component of momentum equation 
orthogonal to it (and to B). 

• Additional drift (diamagnetic) arises in standard FAB form from pressure force. 

• Diagmagnetic drift is opposite for opposite signs of charge (electrons vs. ions). 

Now restore species distinctions and consider electrons plus single ion species i. Quasineu- 
trality says riiqi = —Ueqe- Hence adding solutions 

riegeVe + rij^iVj = {njqi + ng^e)^ -V {pe +Pi) (4-71) 



Hence current density: 

j = -V(pe+P.)A-^ (4.72) 

This is the diamagnetic current. The electric field, E, disappears because of quasineutrality. 
(General case J^jQj^jVj = — V(I]Pj)AB/S^). 

4.6 Reduction of Fluid Approach to the Single Fluid 
Equations 

So far we have been using fluid equations which apply to electrons and ions separately. 
These are called 'Two Fluid' equations because we always have to keep track of both fluids 
separately. 

A further simplification is possible and useful sometimes by combining the electron and ion 
equations together to obtain equations governing the plasma viewed as a 'Single Fluid'. 
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Recall 2-fluid equations: 

UTl ■ 

Continuity (Cj) + ^.{rijVj) = 0. (4.73) 

Momentum (Mj) mjUj + v^-.vj v^- = njqj (E + v^- A B) - Vpj + Fjk (4.74) 

(where we just write F^^ = —Vj^njUij (vj — v^) for short.) 

Now we rearrange these 4 equations (2x2 species) by adding and subtracting appropriately 
to get new equations governing the new variables: 



Mass Density pm 


= rierrie + riinii 


(4.75) 


C of M Velocity V 


= {nerrieVe + riimiVi) / p^ 


(4.76) 


Charge density pg 


= qeUe + qiUi 


(4.77) 


Electric Current Density j 


= qerie^e + qiUi^i 


(4.78) 




= qeiie (ve — Vj) by quasi neutrality 


(4.79) 


Total Pressure p 


= Pe+Pi 


(4.80) 



1st equation: take nie x Ce + rrii x Ci 



(1) ^ + V. {pmV) = Mass Conservation (4.81) 

L/t' 



2nd take q^ x 0^ + qi x Ci 



dpq 



(2) ^ + V.j = Charge Conservation (4.82) 

3rd take Mg + Mj. This is a bit more difficult. RHS becomes: 

n,qj (E + V, A B) - Vpj + = p^E + j A B - V {pe + Pi) (4.83) 
(we use the fact that F^i — Fie so no net friction) . LHS is 

E^i^j- (^+Vi-vji;,- (4.84) 

The difficulty here is that the convective term is non-linear and so does not easily lend itself 
to reexpression in terms of the new variables. But note that since mg « rrii the contribution 
from electron momentum is usually much less than that from ions. So we ignore it in this 
equation. To the same degree of approximation V ~ Vj: the CM velocity is the ion velocity. 
Thus for the LHS of this momentum equation we take 

+ V.-V) V, ^Pmi^^+ V.v") V (4.85) 



75 



so: 



(3) 



dt 



+ V.V V = p,E+jAB- Vp 



Finally we take + to get: 



dt 



+ (v,-V) 



v. = E{^(E + v,AB)--^Vp, 



rrii 



Fjk} 



(4.86) 



(4.87) 



Again several difficulties arise which it is not very profitable to deal with rigorously. Observe 
that the LHS can be written (using quasineutrality niQi + UeQe = 0) as pm§i (^) provided 
wc discard the term in (v.V)v. (Think of this as a linearization of this question.) [The 
(v.V)v convective term is a term which is not satisfactorily dealt with in this approach to 
the single fluid equations.] 

In the R.H.S. we use quasineutrality again to write 



^ rrii 



2 2 



nip 



Heme niTJli 
2 

e 



2 2^i'^i + '^eUe 



E 



nemeUimi 



merrii 



Pm^, (4.88) 



rrii 



1 Ve H Vj} 

merrii q^ 



merrii 



rrii rrie 

-jnemeVe + niVfliWi \ 

(rrii rne\ 

{PmV - — + j} 

\<li Qe J 



{qeneVe + ^i'T'iVj)} 



Also, remembering F^i = -i/ei?T'e"^i(Ve - Vj) = -F^e, 



-Vei ri^qe - rieqi 



rrii 



(Ve - Vj) 



^ _ qe^m^\ . 
Qi rrii J 



So we get 



dt \p„ 



merrii 



PmE + I PmV 



rrii rrie 

, Qi Qe 



j AB 



— Vpe ypi 

rrie rrii 



qeme\_. 

1 h^eij 

Qi rrii I 



(4.89) 



(4.90) 



(4.91) 



Regroup after multiplying by ^^'^ 
E+VAB = 



rrifmi d f j 



qeqi dt \Prn) 
mp rrii 



1 / mj 


me\ . 


— — + 




Pm V 1i 


Qe ) 


\ merrii 


■(- 


) PrnQeqi 





jAB 



qeme 



memi 



qi rrii J qeQiPm 



(4.92) 
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Notice that this is an equation relating the Electric field in the frame moving with the fiuid 
(L.H.S.) to things depending on current j i.e. this is a generahzed form of Ohm's Law. 

One essentially never deals with this full generalized Ohm's law. Make some approximations 
recognizing the physical significance of the various R.H.S. terms. 

Tfl^TtXi ^ I J \ 

— - — -— — arises from electron inertia. 
QeQi at \p„J 

it will be negligible for low enough frequency. 

— I — - -\ ^ I j A B is called the Hall Term. 

Pm \qi qe J 

and arises because current flow in a B-fleld tends to be diverted across the magnetic field. 
It is also often dropped but the justification for doing so is less obvious physically. 

— Vpi term <^ — Vpe for comparable pressures, 
rrii rUe 

and the latter is ~ the Hall term; so ignore qjSJpi/mi. 

Last term in j has a coefficient, ignoring mg/mj c.f. 1 which is 

-. r = — — = T] the resistivity. (4.93) 

qeqi {riimi) qjUe 

Hence dropping electron inertia. Hall term and pressure, the simplified Ohm's law becomes: 

E + V A B = r;j (4.94) 

Final equation needed: state: 

Pen~'^'' +Pinl'^^ — constant. 
Take quasi- neutrality ^ rig (x Ui (x pm- Take ■je = Ji, then 

PP^ — const. (4.95) 

4.6.1 Summary of Single Fluid Equations: M.H.D. 

Mass Conservation : ^ + V {pmV) = (4.96) 

dp 

Charge Conservation : — + V.j = (4.97) 

at 

Momentum : '^"^ + V-V^ V = p^E + j A B - Vp (4.98) 

Ohm's Law : E + V A B = 77 j (4.99) 

Eq. of State : pp^^ = const. (4.100) 
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4.6.2 Heuristic Derivation/Explanation 

Mass Charge: Obvious. 

Mom"" p„||^ + V.vjv= pqE + - Vp (4.101) 

V ^ ^ Electric Magnetic Force PreSSUre 



rate of change of body force on current 

total momentum density 



Ohm's Law 



The electric field 'seen' by a moving (conducting) fiuid is E + V A B = Ey electric field in 
frame in which fluid is at rest. This is equal to 'resistive' electric held r]y. 

Ey = E + VAB = r;j (4.102) 

The pqE term is generally dropped because it is much smaller than the j A B term. To see 
this, take orders of magnitude: 

V.E = p,/eo so Pg^Eeo/L (4.103) 
VAB = //oj l+^^j so aE^jr^B/iioL (4.104) 



Therefore 



PqE eo f B \^ LjiQ I(? I light transit time \ 



jB L y^QaL J B"^ [nQaL'^y \^ resistive skin time J ^ ^ 

This is generally a very small number. For example, even for a small cold plasma, say T^ — l 
eV (o" ^ 2 X 10^ mho/m), L = 1 cm, this ratio is about 10~^. 

Conclusion: the pqE force is much smaller than the j A B force for essentially all practical 
cases. Ignore it. 

Normally, also, one uses MHD only for low frequency phenomena, so the Maxwell displace- 
ment current, dE/c^dt can be ignored. 

Also we shall not need Poisson's equation because that is taken care of by quasi-neutrahty. 



4.6.3 Maxwell's Equations for MHD Use 

V.B = ; VAE=— ; V A B = J . (4.106) 

The MHD equations find their major use in studying macroscopic magnetic confinement 
problems. In Fusion we want somehow to confine the plasma pressure away from the walls 
of the chamber, using the magnetic field. In studying such problems MHD is the major tool. 

On the other hand if we focus on a small section of the plasma as we do when studying 
short-wavelength waves, other techniques: 2-fluid or kinetic are needed. Also, plasma is 
approx. uniform. 
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'Macroscopic' Phenomena MHD 
'Microscopic' Phenomena 2-Fluid/Kinetic 



4.7 MHD Equilibria 

Study of how plasma can be 'held' by magnetic field. Equilibrium =^ V = ^ = 0. So 
equations reduce. Mass and Faraday's law are ~ automatic. We are left with 

(Mom") ^ 'Force Balance' = j A B - Vp (4.107) 

Ampere V A B = // J (4.108) 

Plus V.B = 0, V.j = 0. 

Notice that provided we don't ask questions about Ohm's law. E doesn't come into MHD 
equilibrium. 

These deceptively simple looking equations are the subject of much of Fusion research. The 
hard part is taking into account complicated geometries. 

We can do some useful calculations on simple geometries. 



4.7.1 6'-pinch 




Figure 4.5: ^-pinch configuration. 

So called because plasma currents flow in ^^-direction. 

Use MHD Equations 

Take to be oo length, uniform in z-dir. 

By symmetry B has only z component. 

By symmetry j has only 6- comp. 

By symmetry Vp has only r comp. 
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So we only need 



I.e. 



Force (j A B)^ - (Vp), = (4.109) 

Ampere (V A B)^ = (/xj)^ (4.110) 

i.e. jeB, - |:P = (4.111) 

^ ^ (4.112) 



dr 

lio dr dr 



Eliminate j: - — ^ - |^ = (4.113) 



Solution -^+p = const. (4.115) 




Figure 4.6: Balance of kinetic and magnetic pressure 

-^+p^^^ (4.116) 

[Recall Single Particle Problem] 

Think of these as a pressure equation. Equilibrium says total pressure = const. 

+ p = const. (4.117) 

^-\/^ kinetic pressure 

magnetic pressure 
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Ratio of kinetic to magnetic pressure is plasma 

- ^ (4.118) 

z 

measures 'efficiency' of plasma confinement by B. Want large /3 for fusion but limited by 
instabilities, etc. 

4.7.2 Z-pinch 




Figure 4.7: Z-pinch configuration, 
so called because j flows in z-direction. Again take to be oo length and uniform. 

j=j,e, B^Beeg (4.119) 



Eliminate j: 



or 



Force (j A B)^ - (Vp), = -j,Be - = 



Ampere (V A B)^ - (//j) 

B, d 



ld_ 
r dr 

dp 



[rBe) - Hojz = 



l^oT dr 



(rBe) - ^ = 



IjLoT dr I 2//o 



+ 



p] =0 



(4.120) 
(4.121) 

(4.122) 
(4.123) 



Extra Term Magnetic+Kinetic pressure 

Extra term acts like a magnetic tension force. Arises because B-field lines are curved. 
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Can integrate equation 



bB^ dr 



J a Ho 



+ 



Bl 
2iio 



+ p{r) 



If we choose b to be edge {p{b) — 0) and set a = r we get 




Figure 4.8: Radii of integration limits. 



(4 



2//o 2//o Jr Ho r' 



(4 



Force balance in z-pinch is somewhat more complicated because of the tension force 
can't choose p(r) and j{r) independently; they have to be self consistent. 

Example j — const. 



-— [rBe) = Hojz =^ -Be = r 
r or 2 



Hence 



p(r) 



1 / /^oJ; ' ^ 

2fio V 2 



{b^-r^+ [\r'dr'} 

Jr 



(4 

(4 
(4 



P 




Figure 4.9: Parabolic Pressure Profile. 

Also note Be{b) = ^ so 

-'^66 ^ rT2 2i 
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4.7.3 'Stabilized Z-pinch' 

Also called 'screw pinch', 9 — z pinch or sometimes loosely just 'z-pinch'. 
Z-pinch with some additional as well as Bq 

(Force), - j,Bg - ^ = (4.130) 
d 

Ampere : = l^oje (4.131) 

\l^{rhe)^l^o3. (4.132) 

Eliminate j: 

jio or /lor or Or 

or 

— +^if + ") =» ("34) 

/V^ dr \2no J 

^"e Illy""" (6>+z)+Kinetic pressure 

4.8 Some General Properties of MHD Equilibria 

4.8.1 Pressure &; Tension 

jAB-Vp = : VAB = //J (4.135) 
We can eliminate j in the general case to get 

— (V A B) A B = Vp. (4.136) 

Expand the vector triple product: 

Vp= — (B.V)B- — VS^ (4.137) 
put b = ^ so that VB = WBh = BWh + hWB. Then 

Vp = — {s2(b.V)b + Sb(b.V)S} - -^VS^ (4.138) 
= — (b.V)b- — (V-b(b.V))S2 (4.139) 

= -(b.V)b-vJ— (4.140) 
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Now V_L (^) is the perpendicular (to B) derivative of magnetic pressure and (b.V)b is the 
curvature of the magnetic field line giving tension. 

I (b.V)b I has value ^. R: radius of curvature. 

4.8.2 Magnetic Surfaces 

= B. [j A B - Vp] = -B.Vp (4.141) 

*Pressure is constant on a field line (in MHD situation) . 
(Similarly, = j.p A B - Vp] = j.Vp.) 




p= const, contours 



Figure 4.10: Contours of pressure. 

Consider some arbitrary volume in which Vp ^ 0. That is, some plasma of whatever shape. 
Draw contours (surfaces in 3-d) on which p = const. At any point on such an isoberic surface 
Vp is perp to the surface. But B.Vp = implies that B is also perp to Vp. 




Figure 4.11: B is perpendicular to Vp and so lies in the isobaric surface. 

Hence 

B hes in the surface p = const. 

In equilibrium isobaric surfaces are 'magnetic surfaces'. 

[This argument does not work if p = const, i.e. Vp = 0. Then there need be no magnetic 
surfaces.] 



84 



4.8.3 'Current Surfaces' 



Since j.Vp = in equilibrium the same argument applies to current density. That is 

j lies in the surface p = const. 

Isobaric Surfaces are 'Current Surfaces'. 

Moreover it is clear that 

'Magnetic Surfaces' are 'Current Surfaces'. 

(since both coincide with isobaric surfaces.) 

[It is important to note that the existence of magnetic surfaces is guaranteed only in the 
MHD approximation when Vp 7^ > Taking account of corrections to MHD we may not 
have magnetic surfaces even if Vp 7^ 0.] 



4.8.4 Low P equilibria: Force-Free Plasmas 

In many cases the ratio of kinetic to magnetic pressure is small, /3 << 1 and we can approx- 
imately ignore Vp. Such an equilibrium is called 'force free'. 

j A B = (4.142) 

implies j and B are parallel, 
i.e. 

j = ^(r)B (4.143) 
Current flows along field lines not across. Take divergence: 

= V.j = V.(At(r)B)=/x(r)V.B + (B.V)At (4.144) 
= (B.V)At. (4.145) 

The ratio j/B = fi is constant along field lines. 

II is constant on a magnetic surface. If there are no surfaces, 11 is constant everywhere. 



Example: Force-Free Cylindrical Equil. 



j A B = ^ j = ii{r-)B (4.146) 
VAB = /xj = /x„/x(r)B (4.147) 

This is a somewhat more convenient form because it is linear in B (for specified A*(r)). 

Constant-Ai : V A B = n^jjiB (4.148) 

leads to a Bessel function solution 
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B, = BoJoiHofJ^r) 
Be = BoJiiiiol^r) 



(4.149) 
(4.150) 



for iio/jir > 1st zero of Jo the toroidal field reverses. There are plasma confinement schemes 
with II ~ const. 'Reversed Field Pinch'. 

4.9 Toroidal Equilibrium 

Bend a z-pinch into a torus 




Figure 4.12: Toroidal z-pinch 

Bg fields due to current are stronger at small R side =^ Pressure (Magnetic) Force outwards. 
Have to balance this by applying a vertical field to push plasma back by A B^,. 




Figure 4.13: The field of a toroidal loop is not an MHD equilibrium. Need to add a vertical 
field. 

Bend a ^-pinch into a torus: B(f) is stronger at small R side =^ outward force. 

Cannot be balanced by B^ because no j^. No equilibrium for a toroidally symmetric ^-pinch. 

Underlying Single Particle reason: 

Toroidal ^-pinch has 5^ only. As we have seen before, curvature drifts are uncompensated 
in such a configuration and lead to rapid outward motion. 
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Figure 4.14: Charge-separation giving outward drift is equivalent to the lack of MHD toroidal 
force balance. 



We know how to solve this: Rotational Transform: get some Bq. Easiest way: add j^. Prom 
MHD viewpoint this allows you to push the plasma back by A force. Essentially, this 
is Tokamak. 



4.10 Plasma Dynamics (MHD) 



When we want to analyze non-equilibrium situations we must retain the momentum terms. 
This will give a dynamic problem. Before doing this, though, let us analyse some purely 
Kinematic Effects. 

'Ideal MHD' ^ Set eta = in Ohm's Law. 

A good approximation for high frequencies, i.e. times shorter than resistive decay time. 

E + V A B - 0. Ideal Ohm's Law. (4.151) 

Also 



V A E = Faraday's Law. (4.152) 

Together these two equations imply constraints on how the magnetic field can change with 
time: Eliminate E: 

+ VA(VAB) = +— (4.153) 
This shows that the changes in B are completely determined by the flow, V. 



4.11 Flux Conservation 

Consider an arbitrary closed contour C and spawning surface S in the fluid. 
Flux linked by C is 

* = / B.ds (4.154) 
•J s 

Let C and S move with fluid: 

Total rate of change of $ is given by two terms: 



87 



Contour at time t, t'=t+dt 
Figure 4.15: Motion of contour with fluid gives convective flux derivative term. 



dB 



As 



+ i B.(VAdl) 



c 



Due to changes in B Due to motion of C 

- / VAE.ds- /(VAB).dl 

Js Jc 

- / (E + V A B).dl = by Ideal Ohm's Law. 
J c 



Flux through any surface moving with fluid is conserved. 



(4.155) 

(4.156) 
(4.157) 



4.12 Field Line Motion 

Think of a field line as the intersection of two surfaces both tangential to the field everywhere: 




Field Line 



Figure 4.16: Field line defined by intersection of two flux surfaces tangential to field. 
Let surfaces move with fluid. 

Since all parts of surfaces had zero flux crossing at start, they also have zero after, (by flux 
conserv.). 



Surfaces are tangent after motion 



=^ Their intersection defines a field line after. 

We think of the new field line as the same line as the old one (only moved). 
Thus: 

1. Number of field lines (= flux) through any surface is constant. (Flux Cons.) 

2. A line of fluid that starts as a field line remains one. 



4.13 MHD Stability 

The fact that one can find an MHD equilibrium (e.g. z-pinch) does not guarantee a useful 
confinement scheme because the equil. might be unstable. Ball on hill analogies: 




Stable Unstable Unstable 



Figure 4.17: Potential energy curves 

An equilibrium is unstable if the curvature of the 'Potential energy surface' is downward 
away from equil. That is if ^{M/pot} < 0. 

In MHD the potential energy is Magnetic + Kinetic Pressure (usually mostly magnetic). 

If we can find any type of perturbation which lowers the potential energy then the equil is 
unstable. It will not remain but will rapidly be lost. 

Example Z-pinch 

We know that there is an equilibrium: Is it stable? 
Consider a perturbation thus: 











I 1 

1 V 




J 1 




1 T ^ 



































Figure 4.18: 'Sausage' instability 



Simplify the picture by taking the current all to flow in a skin. We know that the pressure 
is supported by the combination of B^/2fio pressure and ^ tension forces. 
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Figure 4.19: Skin-current, sharp boundary pinch. 
At the place where it pinches in (A) 

Be and ^ increase Mag. pressure & tension increase =^ inward force no longer balance 
by p ^ perturbation grows. 

At place where it bulges out (B) 

Be ^ ^ decrease — Pressure & tension ^ perturbation grows. 

Conclusion a small perturbation induces a force tending to increase itself. Unstable (= 5W < 
0). 

4.14 General Perturbations of Cylindrical Equil. 

Look for things which go like exp[i(A;2; + mO)]. [Fourier (Normal Mode) Analysis]. 




Figure 4.20: Types of kink perturbation. 
Generally Helical in form (like a screw thread). Example: m = 1 7^ z-pinch 

4.15 General Principles Governing Instabilities 

(1) They try not to bend field lines. (Because bending takes energy). Perturbation 
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Figure 4.21: Driving force of a kink. Net force tends to increase perturbation. Unstable. 




■'I'H^^lTt tic 



Figure 4.22: Alignment of perturbation and field line minimizes bending energy. 

(Constant surfaces) lies along magnetic field. 
Example: 6'-pinch type plasma column: 




Figure 4.23: 'Flute' or 'Interchange' modes. 

Preferred Perturbations are 'Flutes' as per Greek columns — > 'Flute Instability' [Better 
name: 'Interchange Instability', arises from idea that plasma and vacuum change places.] 

(2) Occur when a 'heavier' fluid is supported by a 'lighter' (Gravitational analogy). 

Why does water fall out of an inverted glass? Air pressure could sustain it but does not 
because of Ray leigh- Taylor instability. 

Similar for supporting a plasma by mag field. 

(3) Occur when | B \ decreases away from the plasma region. 
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Figure 4.24: Inverted water glass analogy. Rayleigh Taylor instability. 



R2 r2 

2/io 2/io 



(4.158) 



Perturbation Grows. 



(4) Occur when field line curvature is towards the plasma (Equivalent to (3) because of 
V A B = in a vacuum) . 




Figure 4.25: Vertical upward field gradient is unstable. 
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Figure 4.26: Examples of magnetic configurations with good and bad curvature. 
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4.16 Quick and Simple Analysis of Pinches 



^-pinch I I = const, outside pinch 

= No field line curvature. Neutral stability 

2;-pinch V | 5 | away from plasma outside 

= Bad Curvature (Towards plasma) =^ Instability. 

Generally it is difficult to get the curvature to be good everywhere. Often it is sufficient 
to make it good on average on a field line. This is referred to as 'Average Minimum B'. 
Tokamak has this. 

General idea is that if field line is only in bad curvature over part of its length then to perturb 
in that region and not in the good region requires field line bending: 



^ Perturbation prefered by 
-curvature, bends field line. 



Plasma Region 



Good Curvature 



Figure 4.27: Parallel locahzation of perturbation requires bending. 

But bending is not preferred. So this may stabilize. 

Possible way to stabilize configuration with bad curvature: Shear 
Shear of Field Lines 




Figure 4.28: Depiction of field shear. 



Direction of B changes. A perturbation along B at is not along B at Z2 or zi so it would 
have to bend field there — > Stabilizing effect. 

General Principle: Field fine bending is stabihzing. 

Example: Stabilized z-pinch 

Perturbations (e.g. sausage or kink) bend so the tension in Bz acts as a restoring force 
to prevent instability. If wave length very long bending is less. =^ Least stable tends to be 
longest wave length. 

Example: 'Gylindrical Tokamak' 
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Tokamak is in some ways like a periodic cylindrical stabilized pinch. Longest allowable wave 
length = 1 turn round torus the long way, i.e. 



A = 27rR. 



(4.159) 



Express this in terms of a toroidal mode number, n (s.t. perturbation oc expi(n0 + mO): 

4> = jl n = kR. 

Most unstable mode tends to be n = 1. 

[Careful! Tokamak has important toroidal effects and some modes can be localized in the 
bad curvature region {n ^ 1). 



Figure 4.29: Ballooning modes are locahzed in the outboard, bad curvature region. 



/ 



Good 
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